Effects of internal fluctuations on the spreading of Hantavirus 
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Abstract 

We study the spread of Hantavirus over a host population of deer mice using a population dy- 
namics model. We show that taking into account the internal fluctuations in the mouse population 
due to its discrete character strongly alters the behaviour of the system. In addition to the fa- 
miliar transition present in the deterministic model, the inclusion of internal fluctuations leads 
to the emergence of an additional deterministically hidden transition. We determine parameter 
values that lead to maximal propagation of the disease, and discuss some implications for disease 
prevention policies. 

PACS numbers: 87.19.Xx, 87.23. Cc, 05.40.-a 



1 



I. INTRODUCTION 



Hantavirus epidemics have been studied extensively in the biological literature following 
a number of outbreaks in the North American Southwest in the 1990's Q]. The host of this 
infection is the deer mouse, the most numerous mammal in North America. The virus is 
transmitted among deer mice via biting and to humans via contact with their excreta. 

Recently, the Hantavirus has been receiving increasing attention in the physical and 
mathematical literature. A basic population dynamics model was introduced and solved 
by Kenkre P|, the spatiotemporal patterns of the infection were studied by Abramson and 
Kenkre P], Monte Carlo simulations were performed by Aguirre et al. jj], propagating fronts 
of the disease were analyzed by Abramson et al. ^, and the relation between outbreaks of 
the disease and seasonal changes was explored by Buceta et al. jg]. This work has shed light 
on the mechanisms of propagation of the disease among mice, and will hopefully help design 
more effective prevention policies. 

In this paper we go a step further and analyze the effects of the internal fluctuations 
on the propagation of the disease. These fluctuations are inevitable because the mouse 
population is discrete and finite, and they may have profound consequences, as reported by 
Escudero et al. Ql using a generic population dynamics model. 



The basic model introduced in 
and infection. The model reads: 



incorporates birth, death, competition for resources. 



— = (6-c)M--, (la) 

^|i = -eM,-^f£+a(M-A«M,. (lb) 

at K 

where M stands for the total number of mice, M/ for the total number of infected mice, h is 
the birth rate coefficient, c the death rate coefficient, a the infection rate coefficient, and K 
the carrying capacity that characterizes the resources available to the mice and the resulting 
competition. The steady state value of the total mouse population is M = K{h — c). One 
can see that there is a transcritical bifurcation sX K = Kc, with 

[a{b - c)] 

When K < Kc the stable point 

M = K{b-c), Mi = (3) 



has zero infected mice, while when K > the stable point includes a positive number of 
infected mice, 

M = K{b-c), Mi = K{h-c)--. (4) 

(X 

The two rate equations can be thought of as describing two "reactants," M and M/, under- 
going four types of "reactions" with rate coefficients a, 6, c, and respectively. One of 
these conserves the total number of mice (the infection), while the other three (birth, death, 
competition) do not. Note that Eq. (jlap depends only on the latter three, whose effect on 
the total population can thus be studied separately from the issue of infection. Infected 
pregnant mice produce Hanta antibodies that keep their foetus free from the infection. Con- 
sequently, there is no birth term in Eq. ()lb|) . Also, there is no recovery term in the model 
because mice become chronically infected with the virus. 

The analysis of the internal fluctuations in the mouse population due to the discrete 
and finite sizes of the populations requires the generalization of the mean field model to a 
stochastic description, e.g., a master equation. In this paper we start with such a master 
equation, but at the very outset we outline some reasonable approximations that lead to a 
mathematically tractable model. 

A full master equation description of the problem would involve P(n, n/,t), the prob- 
ability distribution function for there to be n total mice and nj infected mice at time t. 
We find this full master equation to be analytically intractable. We therefore break the 
problem up into two parts as follows. First, we formulate a master equation for the reduced 
probability distribution function P{n,t) associated only with the mean field equation ()la|). 
This master equation (which is not infiuenced by the infection) is tractable, as we shall 
see. We then argue that the fiuctuations in the infected mouse population arise from two 
sources. One is the dependence on M in Eq. (jlbj) and the fact that this total population 
fiuctuates. Having solved the master equation associated with M, we are able to incorporate 
these fiuctuations into the stochastic description of infected mice. We will show that the 
effects of these fiuctuations may be profound, especially when the mean mouse population 
is not too large, and may lead to unexpected consequences. These are the new features that 
we are particularly interested in exploring. The other arises from the additional inherent 
fiuctuations in the number of infected mice due to the fact that this population is also finite 
and discrete. These are especially important when the number of infected mice is small, but 
we do not include them explicitly in our equations, again because of tractabihty problems. 
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This is not as serious as one might think because we do know their consequences, which can 
also be profound (as we have shown in [3]): these fluctuations may cause a small popula- 
tion of infected mice to disappear entirely. In other words, if one is in a regime where the 
population of infected mice is small in the absence of these fluctuations, consideration of 
these fluctuations might eliminate this population entirely. Thus, results obtained without 
consideration of these fluctuations can be thought of as an upper bound on the number of 
infected mice. At worst one would be overestimating the presence of infected mice in the 
regime where the number of infected mice is in any case small or zero. 

In Sec. nil we present the stochastic model for the total mouse population. Section IIIII 
deals with the stochastic model for the infected mouse population, and in Sec. [IVlwe discuss 
the results of the analysis. We summarize our conclusions in Sec. [V] 



II. STOCHASTIC MODEL FOR TOTAL MOUSE POPULATION 

The master equation for the total mouse population is easily written down if we think 
explicitly of the "reactions" contributing to Eq. ()la|). They are births. 



deaths, 



and competition for resources. 



M M + M, (5) 



M A 0, (6) 



M + M^-^ M. (7) 



The master equation describing these processes is 
dP{n,t) 



dt 



b [{n - l)P{n - 1, t) - nP{n, t)] +c[{n + l)P{n + 1, t) - nP{n, t)] 
+K-^ [(n + l)nP{n + 1, t) - n{n - l)P{n, t)] . (8) 



^ as it stands, but it is amenable to a system size expansion as 
0]. A system size expansion is appropriate when the system 



This equation is not tractab^ 
introduced by van Kampen 
is "large" or, as in our case, the species under consideration numerous. It is important to 
note that this implies that the steady state solution for M given in Eqs. (jH]) and (jH), which 
we expect the mean of the stochastic solution to reproduce, must therefore be "large", that 
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is, K must be proportional to the system size. The ratio 



(9) 



where VL is the size of the system, must be essentially independent of the system size for this 
analysis to be appropriate. To implement a system size expansion we thus write the third 
coefficient on the right of Eq. (jHJ as = d/Q. 

The implemention of a system size expansion requires several steps. First, although n 
is a discrete variable, we can represent the discrete changes in n via an infinite series of 
derivatives in which n is treated as a continuous variable: 



/ j=0 



[±iy 



fin). 



(10) 



This exact relation allows us to rewrite the master equation (jH)) as 
dP{n,t) 



dt 



1)^ ^ 1 

=1 •' 



j ! dn^ 



nP(n,t) + ^^i|^n(n-l)P(n,t). (11) 



=1 



Next, one makes the heuristic assumption that one can perform the change of variables 

n n(j){t) + n^^^z, (12) 

where 0(t) is the mean value of the mouse population density and z represents the fiuctua- 
tions around the mean. We then define the probability distribution 

P{n,t) 



Applying the standard chain rule 

'dP{n,ty 



p{z,t) 



dpiz,t) 
dt 



dz 



together with the relation which follows from the change of variables relation, 

d(t){t) 



fdz^ 

\dt 



dt 



(13) 



(14) 



(15) 



we can (after multiplying through by Q}/'^) rewrite the master equation in terms of the new 
distribution, 



dp{z,t) 
dt 



d(f) dp{z, t) 
dt dz 



^ ]\ dzi ^ 



1 

j! dz^ 



(fi0 + l]i/2^)p(z,t) 



j! dzi 



(16) 



This equation is still exact. 

In the large system size limit there are three divergent terms in Eq. ()16|) proportional to 
Vt^^'^dp/dz that must cancel, that is, we must require that 

^ = (6 - c)</> - dct.\ (17) 

Note that this exactly corresponds to Eq. (fTaj) . that is, (f){t) is indeed the mean population 
density. In the steady state we thus have 

Using this result in the surviving terms in Eq. ()16|) in the large VL limit then leads, in the 
steady state, to the equation 

dt ^ ^ dz ^ d dz'' 

^d{zpl^b^_ (19) 
dz d dz'^ 

This is the steady state limit of a Fokker-Planck equation for the probability density of the 
stochastic variable z. Its solution with natural boundary conditions at ±oo is given by 



a Gaussian distribution centered at zero and of width proportional to yhjd. This in turn 
implies that the total number n of mice also has a Gaussian distribution whose mean is the 
mean number of mice predicted by the deterministic model and whose width is proportional 
to ^/Kh. 

The internal fluctuations thus do not alter the behavior of the total number of mice in any 
dramatic way. They simply lead to a Gaussian distribution around the deterministic mean 
whose width increases with increasing birth rate and increasing carrying capacity. However, 
as we will see in the following sections, the consequences of this distribution on the number 
of infected mice can be unexpected. 



III. STOCHASTIC MODEL FOR INFECTED MOUSE POPULATION 

We now return to the infected mouse population, whose evolution is described in mean 
field by Eq. ()lb|) . While we are ignoring the internal fluctuations that arise from the fact 

6 



that this population is finite and discrete (as discussed earher), we do wish to provide a 
stochastic description that incorporates the effects of the fluctuations in the total mouse 
population n(t). Our results of the previous section indicate that in the steady state we can 
think of n{t) as a stochastic variable, 

n{t)=K{b-c) + 6n{t) = M + 6n{t), (21) 

where the fluctuations Sn{t) have zero mean and are generated from the Ornstein-Uhlenbeck 
stochastic differential equation 

-{b - c)6n + ^/2Kb{b - c)C{t). (22) 



dt 



Here ^{t) is zero-centered 5-correlated Gaussian noise of unit intensity, {^{t)^{t')) = 5{t — t'). 
In the stationary state the correlation function of the fluctuations is then 

{5n{t)Sn{t')) = i^6e-(^-^)l*-*'l. (23) 

We include these fluctuations in Eq. ()lb|) by replacing n with K{b — c) + 6n. The resulting 
stochastic differential equation reads 

^ = [aK{b -c)-b]nj- an\ + ^J_^' n,C(t), (24) 

where C^it) is an Orstein-Uhlenbeck process with zero mean and correlation function 

(C(t)C(t')) = ^K&-c)e-(^-^)l*-*'l (25) 

This is our basic stochastic equation for the infected population. The intensity of these 
fluctuations is determined by the width of the total mouse population distribution. The 
correlation time = {b — c)~^ is a measure of the time it takes a total population diminished 
by fluctuations to recover. 

Because ({t) is a "colored noise" with flnite correlation time, the exact solution of the 
problem is not known. In particular, there is no exact Fokker-Planck equation 

for the probability distribution function P{nj,t) that the number of infected mice is nj at 
time t. A number of approximate Fokker-Planck equation schemes can be found in the 
literature some of which have the virtue of becoming exact in both the limits Tc ^ 
and Tc — > cxD. Since these theories lead to a qualitatively similar panorama of possibilities. 
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The resulting effective 



we apply the simplest of these theories, developed by Fox 
Fokker-Planck equation is 



where 



G{nj) = [aK{b — c) — b]nj — anj, 



(26) 



(27) 



and 



, , a-K-^ 

g{ni) = ni, 

y/b — C 

Kb{b - c) 



D{nj) 



b — c + aril 



The stationary solution of this equation is 



P{ni) = N [1 + -^Jii ) rii 



{b - c)K[a{b - c)K - b] 
b{aK-lf 



(28) 
(29) 



X exp 



Ka^n] ^aK[aK{b- c) -2b + c]ni 



(30) 



2b{aK-iy b{aK-iy 
where N is the normalization constant. In the next section we analyze and comment on the 
interesting features of this solution. 



IV. RESULTS FOR INFECTED MOUSE POPULATION 

The first point to note is that the mean number of mice, 

POO 

Ml = {nj) = / dnjniP{nj), (31) 







is exactly as predicted in the mean field theory, cf. Eqs. However, here there are 

distributions of infected mouse populations underlying this mean, and our interest lies in the 
different shapes of these distributions in different parameter regimes, and in the additional 
information beyond the mean contained in these distributions. The distributions, described 
in detail below, are sketched in Fig.Q where we present the phase diagrams of the system in 
{K, a) space for fixed b and c. The values of b and c have been chosen to match those used 



in earlier Monte Carlo simulations 
that are observed in the system. 



and/or to make most evident the different behaviors 
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FIG. 1: Phase diagrams for the infected mouse population in {K,a) space, described in detail in 
the text. 

Various phase boundary lines are shown in Fig. ^ The sohd curves in both panels 
are the curves K = K^. When K < Kc the probability distribution (jHUj) can not be 
normalized because it has a nonintegrable singularity at n/ = 0. Since this is a fixed point 
of the d;piamics, the probability distribution must be interpreted as a 5 function centered 
at zero The insets at the lower left of each panel are a schematic of this behavior. 
When K crosses the Kc curve there is still a divergence at n/ = 0, but the probability 
distribution becomes integrable and hence normalizable. The most probable value for the 
infected population is still zero, but nonzero values now have a finite probability and the 
mean value of infected mice is positive. The lower right hand insets in both panels are 
sketches of this behavior, which persists as a function of increasing K until the carrying 
capacity reaches a second critical value, 

2h 

j^* ("321 

' ~ h[2a - (6 - c)] + ^h{h-c)[h^ -c{Aa + h)] 

When K crosses the curve K = K* the divergence in the distribution (jSOI) disappears 
and the most probable number of infected mice moves to finite values. The curve has a 
divergent asymptote at a = 0, but its behavior as a function of a otherwise depends on 
the other parameters. If 6 < 2c then K* also diverges at a = Oc = (6 — c)'^/b. This is the 
situation in panel (a) of Fig. ^ On the other hand, if 6 > 2c then K* is complex when 
a > ttc = h{h — c)/4c, thus producing the abrupt vertical boundary seen in panel (b). This 
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FIG. 2: Ratio of the dispersion to the mean number of infected mice for b = 0.5 and c = 0.2. 

second case corresponds to the parameters b and c in the Monte Carlo simulations of Aguirre 
et al. [4j. In either case, within the region enclosed by the curve (dashed curve in the 
figure) the probability distribution goes to zero at the origin and has a maximum at a finite 
value of n/, as shown in the upper left sketches in both panels. Both the average number of 
infected mice and the most probable number of infected mice are now positive. Note that we 
have labeled the rightmost value of a on the dashed curve as Oc, whether it is an asymptote 
as in panel (a) or the abrupt ending point of the curve as in panel (b). 

There is an additional transition curve, more subtle than the other two, drawn as the 
dotted curves in both panels in Fig. Q We denote this transition curve as K** . On the 
upper right hand in both panels is the sketch of the probability distribution in this region. 
Here the probability distribution diverges at zero, but another maximum develops at a finite 
number of infected mice. This maximum is found as the finite positive root of the derivative 
condition dP{ni) / drij = 0. The dotted curves in the phase diagrams indicate the location 
of this transition. 

A number of points about these results deserve special highlighting. The particular be- 
havior just described for large K and a (divergence at the origin and also another maximum) 
is entirely due to the fact that the internal fluctuations are colored. The correlation time of 
these fluctuations is (6 — c)^^, and the color has arisen naturally and not as an additional 
assumption. It is interesting to note that the Monte Carlo simulation results of Aguirre et 
al. 0] exhibit a number of features that might be related to the results that we have derived 
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here. One is that in their simulations the number of infected mice as a function of K jumps 
discontinuously from zero to a finite number (whereas the mean field value does not). They 
note that an explanation for this result lies in the fact that the number of mice is discrete 
and that one whole mouse is needed to obtain a result different from zero (although their 
jump is much greater than unity). In our continuous language, the behavior they observe 
might reflect the abrupt transition between the 5-function distribution (or the one with a 
maximum at the origin) to the one with a zero probability density of no infected mice as K 
increases. Their simulations use the value a = 10^^. To support this argument further, we 
have plotted in Fig.|21the ratio of the dispersion to the mean for these parameter values. The 
dispersion is of the order of the mean and, near the transition value Kc, the ratio actually 
diverges. We have elsewhere pursued the argument that a possible criterion for the likely 
extinction of a species is precisely that the dispersion be of the same size as the mean [7|. 
The substantial width of the distribution might make itself apparent in a simulation through 
the high likelihood of absence of the infected species. For comparison, we have also plotted 
the ratio of the dispersion to the mean for a value of a > Oc. The fluctuations are now 
decidedly smaller, even though we are in a regime of far fewer infected mice on average (as 
indicated by the values of K). It would be interesting to see whether the size of the jump in 
the infected mouse population at the transition would decrease in a Monte Carlo simulation 
with a > ttc- 

Finally, we comment on three last points. One concerns the validity of the system size 
expansion. The total mean number of mice in the population is K{b — c), and the system size 
expansion is valid if this number is in some sense sufficiently large (the expansion is valid if 
the neglected terms are small). While we have not explicitly checked the validity, in most of 
the regimes under discussion the number of infected mice is an order of magnitude greater 
than unity. Our second point is to stress that the fluctuations that lead to the distributions of 
infected mice are entirely due to the discrete and finite character of the total number of mice. 
And yet, while the ratio of the width of the distribution of total mice to the mean number 
of total mice, au/M, is small in most of the phase diagram {au/M = {b/K)/{b — c)), 
the width of the distribution induced in the number of infected mice is relatively large (of 
0{> 1)) in most of the diagram. The third point is a reminder that this theory has not 
included the fluctuations caused directly by the fact that the number of infected mice is 
discrete and finite. These fluctuations would further broaden the distributions. 
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V. CONCLUSIONS 



We have considered the effect of the internal fluctuations in the total mouse population 
on the number of infected mice. Although these fluctuations cause no dramatic effects in the 
total mouse population, yielding a Gaussian distribution of relatively small variance for a 
sufficiently large population, these fluctuations have a rather strong effect on the distribution 
of infected mice. Because the fluctuations are not "direct" but instead appear indirectly 
through the coupling between infected and uninfected mice, the fluctuations necessarily 
and naturally appear as colored in the equations that describe the evolution of infected 
mice. This leads to a variety of effects beyond those that would be caused by simple white 



noise 



IG 



The mean infected population in this model is exactly that predicted in mean 



field. However, while the mean field model predicts one critical value of the carrying capacity 
parameter (Kc) such that below this value there is no infection and above this value there 
is, the stochastic model leads to three critical values [K^, K* and K^*). The first, which 
occurs at the same critical value as that of the deterministic model, here corresponds to a 
transition between a state with no infected mice to an intermediate state in which the most 
probable state is still one with no infected mice but with a finite probability of infection. 
The second describes a transition between this intermediate state and the outbreak state, 
where the probability distribution that there is no infection goes to zero. The intermediate 
state {Kf. < K < K*) displays different behaviors depending on the parameter values. In 
particular, in some parameter ranges the intermediate state has very few infected mice. 
We argued that the inclusion of the internal fluctuations in the infected mouse population 
(which was not considered due to analytic difficulties) would probably lead to extinction of 
this small number of infected mice. This then means that the effective transition between 
non-epidemic and epidemic states may occur at K* rather than at K^. We also identified 
another transition curve, K**, beyond which the probability diverges at zero but where 
another maximum develops at a finite number of infected mice. 

These features may be useful in the design of more effective prevention policies. For 
instance, an increase in the effective annihilation rate of the mice (by either increasing the 
death rate or decreasing the birth rate or both) might help because it increases the relative 
size of the region in parameter space in which the infected mouse population distribution 
has a divergence at the origin (the state of no infected mice). The most effective strategy 
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for the control of Hantavirus outbreaks is the reduction of the carrying capacity X so as to 
cross from one regime to another with a higher probabihty of no infected mice. 
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